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1  Introduction 


During  the  contract  period,  we  conducted  a  study  of  the  feasibility  of  controlling  the  input  impedance  of  a 
microstrip  element  by  dynamically  changing  its  loading.  The  motivation  for  this  work  is  the  possibility  of 
using  such  elements  in  a  scanned  array.  By  changing  the  loading  of  the  elements  appropriately,  one  could 
alter  the  active  array  impedance  of  the  elements  perhaps  compensating  to  some  degree  for  the  onset  of  scan 
blindness.  While  this  report  did  not  establish  or  refute  the  feasibility  of  the  ultimate  application,  it  definitely 
establishes  the  feasibility  of  controlling  input  impedance  of  a  single  element  by  using  PIN  shorting  diodes. 

The  basic  theory  of  loaded  elements  is  briefly  reviewed  in  the  next  section.  Those  features  of  the  theory 
pertinent  to  the  current  application  are  discussed  in  more  detail  in  the  third  section.  The  theoretical  results 
and  the  corresponding  experimental  results  are  also  presented  in  this  section.  Conclusions  and  suggested 
future  research  are  given  in  the  last  section.  Additional  details  of  the  analytical  methods  used  can  be  found  in 
Appendix  A.  Many  of  these  theoretical  results  were  just  recently  developed  and  are  not  presented  elsewhere 
in  the  open  literature.  In  particular,  the  extraction  of  quasistatic  terms,  detailed  in  Appendix  A,  makes  the 
interactive  design  of  these  types  of  loaded  elements  feasible.  This  is  because  the  determination  of  the  resonant 
frequency  of  a  loaded  element  is  a  non-linear  problem  requiring  the  iterative  solution  of  the  characteristic 
equation.  Efficiency  in  computation,  therefore,  is  critical — and  the  quasistatic  analysis  presented  is  highly 
efficient. 

2  Theory  of  Loaded  Microstrip  Antennas 

A  brief  overview  of  the  theory  of  loaded  microstrip  antennas  is  presented  in  this  section.  A  more  detailed 
discussion  can  be  found  in  [1]. 

For  design  purposes,  one  needs  a  simple  theory  which  sheds  light  on  the  physical  mechanisms  of  the 
microstrip  antenna.  Even  if  the  theory  is  not  the  most  quantitatively  accurate,  it  is  still  vital  to  understand 
its  physics  if  one  is  going  to  do  creative  things  with  the  element.  One  theory  that  has  proven  to  be  physically 
illuminating,  and  relatively  accurate  as  well,  is  the  cavity-model  analysis.  It  is  this  analysis  that  is  used  as 
a  basis  for  design  of  loaded  microstrip  elements. 

A  typical  microstrip  element  is  illustrated  in  Fig.  1.  The  fundamental  assumption  made  in  cavity-model 
analysis  is  that  the  distribution  of  the  field  between  the  patch  and  the  ground  plane  is  approximately  the 
same  as  that  resulting  in  a  cavity  formed  by  closing  the  edge  of  the  patch  to  the  ground  plane  as  illustrated 
in  Fig.  2.  The  amplitude  of  the  field  distribution  in  the  antenna  and  in  the  cavity  are  different.  The  method 
of  determining  this  amplitude  will  be  discussed  later.  For  the  purpose  of  determining  the  resonant  frequency 
of  the  loaded  element,  the  antenna  will  be  assumed  to  be  the  same  as  the  corresponding  cavity. 


2.1  Unloaded  microstrip  cavity 

In  order  to  analyze  the  loaded  cavity  ultimately  of  interest  in  this  report,  one  must  first  understand  the 
analysis  of  an  unloaded  cavity.  Suppose  one  impresses  a  uniform,  unit,  vertical,  filamentary  electric  current 
as  illustrated  in  Fig.  2.  The  resulting  cavity  electric  field  will  be  vertically  directed.  Let  the  vertically  directed 
electric  field  produced  be  denoted  by  G(r|r'),  where  r  is  the  two-dimensional  position  vector  locating  the 


observation  point  within  the  cavity,  and  r'  is  the  corresponding  position  vector  locating  the  source  current. 
Since  this  field  is  produced  by  a  concentrated  unit  source,  the  field  G(r|r')  can  be  thought  of  as  a  Green’s 
function.  The  boundary-value  problem  that  this  Green’s  function  satisfies  is 

(VJ  +  k2)G(r\r')  =  jup6(r  -  r'),  r  €  S, 

VG(r|r')  •  n  =  0,  r  6  55, 

where  5  denotes  the  set  of  points  on  the  patch,  55  denotes  the  boundary  of  5,  and  n  is  the  outward-pointing 
unit  normal  on  55.  The  solution  of  this  boundary-value  problem  is 


,  /\  ^  Wmnyr  Vm, 

G(r|r  )  =  jw/i  -p 


where  A  is  the  set  of  all  pairs  of  non-negative  integers,  (m,n),  and  ^mn(r)  is  the  mnth  mode  of  the  cavity 
satisfying 

VJV>mn(r)  =  -k^ipmnir),  r  6  S, 

V^mn(r)  •  n  =  0,  i"  €  55. 

The  series  for  G(r  |r')  diverges  as  r  — *  r'.  This  is  because  the  inductive  reactance  associated  with  a 
filament  is  infinite.  In  order  to  obtain  the  input  reactance  of  an  actual  feed,  one  must  integrate  this  Green’s 
function  over  the  actual  feed  current  distribution.  This  integration  is  typically  a  surface  integral  over  the 
outer  periphery  of  the  feed  structure.  For  the  purposes  of  this  report,  it  is  assumed  that  all  source  and  load 
currents  are  distributed  uniformly  on  a  circular  cylinder  of  diameter  d.  The  “location”  of  the  source  or  feed 
is  the  center  of  the  cylinder’s  cross  section.  The  input  impedance  of  the  element  can  be  found  by  placing 
the  observation  point  at  the  center  of  the  cross  section  of  the  feed,  and  by  integrating  the  Green’s  function 
over  the  surface  current.  When  this  is  done,  one  obtains  an  input  reactance  of 

XM  =  Xr  +  Xr, 


Xr  =  —jupt 


jfc2  _  jfc2 
*  *mn 


is  the  resonant  reactance  and  X(  is  the  feed  reactance.  The  parameter  t  is  the  thickness  of  the  dielectric 
under  the  patch.  The  pair  (m,  n)  represents  the  indices  of  the  “resonant  mode.”  For  the  case  of  loaded 
microstrip  elements,  the  resonant  mode  selected  for  determining  XT  is  the  mode  of  the  unloaded  cavity  that 
resonates  at  a  frequency  closest  to  that  of  the  mode  of  interest  in  the  loaded  cavity.  The  feed  reactance 
is  the  contribution  to  the  total  reactance  from  all  of  the  non-resonant  modes,  principally,  the  higher-order 
modes  of  the  element.  The  feed  reactance  is  typically  inductive  and  varies  slowly  with  frequency.  Of  course, 
if  the  frequency  varies  too  far  from  the  resonant  frequency  associated  with  the  indices  (m,  n),  the  feed 
reactance  will  contain  resonant  components  as  well.  The  circuit  model  associated  with  this  is  illustrated  in 
Fig.  3.  The  resonant  mode  of  the  loaded  element  can  be  conveniently  expressed  in  terms  of  the  modes  of  the 
corresponding  unloaded  cavity  as  described  below. 


Xi 


Figure  3:  Circuit  model  for  the  input  impedance  Zm  of  a  microstrip  element  in  the  vicinity  of  one,  well 
isolated  resonant  frequency. 


2.2  Loaded  microstrip  cavity 

The  resonant  frequency  of  a  loaded  cavity  is  that  frequency  at  which  the  cavity  can  support  a  source- 
free  cavity  mode.  At  that  frequency,  in  general,  a  non-sero  current  will  flow  through  the  multiple  loads  of 
the  cavity  even  though  no  source  current  is  applied.  A  non-eeparable  loaded  cavity  can  be  converted  to 
a  separable  cavity  by  using  the  equivalence  principle.  If  the  loads  are  removed  and  equivalent  electric  and 
magnetic  currents  are  impressed  on  the  surface  bounding  the  now  absent  loads,  then  the  field  produced  within 
the  now  unloaded  cavity  is  the  same  as  that  of  the  original  loaded  cavity.  It  turns  out  that  the  magnetic 
currents  called  for  by  the  equivalence  principle  are  typically  negligible  |l).  The  total  electric  current  impressed 
at  the  location  of  each  load  is  the  total  source-free  load  current  that  Sows  through  the  load.  The  unloaded 
cavity,  in  contrast  to  the  loaded  cavity  is  now  driven  by  the  load  currents,  although  the  load  currents  are, 
as  yet,  undetermined.  To  understand  how  to  determine  the  load  currents  and  to  better  understand  the  role 
played  by  load  location,  consider  first  an  element  loaded  by  just  a  single  load. 

2.2.1  Single  load 

Suppose  that  the  load  has  a  reactance  of  X\,.  Then  the  network  model  of  the  undriven  loaded  cavity 
is  that  illustrated  in  Fig.  4.  A  source-free  mode  is  desired.  That  is,  one  seeks  to  produce  a  non-sero  load 
current  with  no  source.  This  is  only  possible  if  the  network  is  resonant.  Summing  the  voltage  around  the 
closed  loop  yields  the  condition  for  resonance, 

Xr  +  X,  +  XL  =  0, 


or 


X,  =*  -X,  -Xu. 

To  best  interpret  this  characteristic  equation,  consider  the  simple  and  very  practical  case  of  a  short- 
circuit  load.  In  this  case,  Xi  «  0.  Figure  5(a)  contains  a  plot  of  Xt  vertut  frequency  for  a  rectangular 
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Figure  4:  Circuit  model  of  &  cavity  loaded  with  a  tingle  load  reactance  of  Xt. 


microstrip  antenna  for  three  different  load  locations  as  illustrated  in  Fig.  5(b).  It  also  contains  plots  of  the 
negative  of  the  feed  reactance  corresponding  to  each  feed  location.  The  intersection  of  corresponding  curves 
is  the  resonant  frequency  of  the  shorted  microstrip  element. 

The  resonant  mode  in  Fig.  5  is  the  (0, 1)  mode,  the  mode  most  commonly  used  in  microstrip  elements. 
For  a  rectangular  element  of  dimensions  a  x  6,  the  normalized  voltage  distribution  at  observation  point  (z,  y) 
for  the  (0,1)  mode  is 

v(x,y)  =  cos  (y)  . 

This  is  zero  on  the  center  line,  y  =  6/2,  of  the  patch.  Such  a  line  of  zero  field  is  called  a  nodal  line.  A  point 
at  which  the  voltage  is  stationary  is  called  an  anti-nodal  point.  All  points  on  the  long  edges  of  the  patch 
illustrated  in  Fig.  5  are  anti-nodal  points.  When  the  load  current  is  impressed  near  a  nodal  line,  the  resonant 
reactance  curve  becomes  very  skinny  as  illustrated  by  case  3  in  Fig.  5(a).  In  contrast,  when  the  load  current 
is  impressed  near  an  anti-nodal  point  such  as  in  case  1  of  Fig.  5(b),  then  the  XT  curve  becomes  fat.  Since 
the  voltage  is  stationary  near  an  anti-nodal  point,  the  Xr  curve  does  not  change  much  as  one  moves  the  load 
in  the  vicinity  of  the  anti-nodal  point.  Thus  the  curve  illustrated  in  case  2  of  Fig.  5(a)  is  not  much  different 
from  the  case  1  curve. 

This  variation  of  the  ,Yr  curve  would  predict  that  the  resonant  frequency  should  monotonically  increase 
as  the  short  is  moved  from  the  nodal  line  to  an  anti-nodal  point  if  the  feed  reactance  were  to  remain  constant. 
However,  as  shown  in  [2],  the  feed  reactance  is  anything  but  constant  with  the  feed  position.  In  fact,  the 
feed  reactance  increases  very  rapidly  as  the  feed  point  approaches  the  edge  of  the  patch.  If  the  edge  is  an 
anti-node,  then  the  Xr  curve  does  not  change  much  as  the  load  is  moved  toward  the  edge,  but  the  —X(  curve 
is  moving  rapidly  downward  in  the  graph  in  Fig.  5  so  that  the  resonant  frequency  of  the  shorted  element 
decreases  for  loads  near  the  edge.  This  makes  sense  because  as  the  feed  reactance  increases,  it  becomes 
more  like  an  open  circuit  and  thus  its  loading  effect  is  reduced.  A  typical  theoretical  frequency  versus  short 
position  variation  is  illustrated  in  Fig.  6.  The  actual  measured  variation  of  resonant  frequency  with  the 
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Figure  5:  (a)  Illustration  of  Xr  and  —Xt  fenu  frequency  for  the  three  feed 
(b)  Rectangular  patch  showing  three  different  feed  locations. 


location  of  the  short  agrees  well  with  the  predicted  results.  The  locat  on  of  the  shorting  pin  is  indicated  by 
the  normalized  shorting  pin  location  defined  as 


If  one  uses  other  types  of  load  reactances,  the  resonant  frequency  of  the  loaded  element  will  be  modified 
although  the  effect  of  the  feed  reactance  is  always  present .  Examples  of  other  applications  of  loaded  microstrip 
elements  can  be  found  in  [3],  [4].  For  the  current  application  of  controlling  input  impedance,  one  should  not 
use  just  a  single  load.  The  reason  why  will  be  explained  in  the  next  section. 

2.2.2  Multiple  loads 

For  multiple  loads,  one  cannot  use  the  simple  graphical  picture  presented  above  for  predicting  the  behavior 
of  a  multiple-loaded  element  although  the  principles  learned  from  this  picture  still  apply.1  For  multiple  loads, 
the  characteristic  equation  is 


Z  11  +  zil) 

Z\ 2 

Z\N 

Z 12 

z22  +  z[7)  . 

Z2N 

=  0, 

(2) 

Z\N 

Z-IN 

..  ZNN  +  ziN) 

where  Zij  is  the  z-parameter  between  the  t*h  and  j'th  port  of  a  multiple-port  microstrip  cavity,  and  is 
the  load  impedance  connected  to  the  Ith  port.  The  z-parameters  can  be  computed  using  equation  (1)  as 
modified  by  integrating  over  the  feed  current  distribution.3  Once  the  resonant  frequency  is  determined,  the 
next  step  is  to  compute  the  modal  field  distribution  of  the  loaded-cavity  resonant  mode. 

To  determine  the  field  distribution  of  the  resonant  mode,  one  must  first  assume  a  unit  (or  any  other 
non-zero  current)  at  one  of  the  ports  and  then  compute  the  resulting  currents  at  the  remaining  loads.  Let 
Z.  be  a  submatrix  excluding  the  first  row  and  column  of  the  matrix  in  equation  (2).  Let  Zi,  be  the  first 
column  of  the  matrix  in  equation  (2),  excluding  the  first  row.  Let  I  be  the  column  matrix  of  port  currents 
excluding  the  current  in  port  1, 

'  h  ' 
h 

1  = 

.  . 

Then,  if  the  load  current  in  port  1  is  I\  =  1,  the  currents  in  the  remaining  ports  are 

i  =  -z.-%„ 

1  Actually,  one  could,  in  principle,  apply  the  graphical  analysis  sequentially.  That  is,  first  find  the  effect  of  a  single  load. 
Then  create  a  new  graph  in  which  Xt  is  the  resonant  reactance  and  Xf  is  the  feed  reactance  of  the  linfly-loaied  pttck.  From 
these,  the  resonant  frequency  of  the  doubly-loaded  element  could  be  constructed. 

3 There  is  more  to  the  computation  than  what  has  been  outlined  here.  A  very  efficient  analysis  has  been  devised  for  computing 
the  z-parameters  in  which  the  series  is  highly  accelerated.  The  details  of  this  are  included  in  Appendix  A. 


(3) 


The  resulting  field  distribution  is 

a  r=i  mn 

where  ^m»(r)  and  ifcmn  are  the  modal  fields  and  corresponding  resonant  wavenumbers  in  the  unloaded  cavity. 
Strictly  speaking,  this  expression  is  only  valid  for  observation  points  r  outside  of  the  load  regions.  That  is, 
if  the  diameter  of  the  Ith  load  is  <U,  then  |r  —  r  j|  must  be  larger  than  di/2. 


The  input  impedance  of  the  loaded  element  is  computed  by  analyzing  the  cavity  with  the  actual  dielectric 
replaced  with  an  artificially  lossy  dielectric — one  whose  losses  account  for  the  radiative  and  ohmic  losses  of 
the  original  microstrip  antenna.  The  effective  loss  tangent  of  this  dielectric  is  Sea  =  1  /Q  where  Q  is  the 
quality  factor  of  the  antenna.  It  is  given  by 

f  i  1  1  IN"1 

Q~  +  +  <3n, +  Q<J  ’ 

where  Qnd  is  the  radiative  quality  factor,  Qtw  is  the  quality  factor  associated  with  the  surface  wave,  Qm  is 
the  quality  factor  associated  with  metal  heating  loss,  and  Qd  is  the  dielectric  quality  factor.  These  are  given 
by 


Qnd  = 


2 uWE 
Pn a  ’ 


Q,  w  = 


2 u,WE 


where  Pnd  is  the  time-averaged  power  radiated  by  the  antenna,  P„  is  the  time-averaged  power  carried  away 
by  the  surface  wave,  A  is  the  skin  depth  in  the  metal  cladding,  6  is  the  loss  tangent  in  the  actual  dielectric, 
and  We  is  the  average  energy  stored  in  the  electric  field  at  the  resonant  frequency  of  the  cavity 


The  parameters  in  the  expressions  for  the  quality  factors  would  appear  to  depend  on  the  field  under  the 
patch.  Actually,  the  quality  factors  depend  only  on  the  distribution  of  the  mode  for  which  the  quality  factors 
are  being  computed.  The  field  could  be  multiplied  by  any  arbitrary  constant  without  affecting  the  quality 
factor.  At  the  resonant  frequency  of  the  ♦(r)  of  the  loaded  element,  the  field  is  strongly  dominated  by 
4(r)  itself.  Thus,  one  can  obtain  an  approximation  for  the  quality  factors  by  assuming  the  field  distribution 
below  the  patch  is  'k(r),  with  the  contributions  of  the  other  non-resonant  modes  of  the  loaded  patch  ignored  . 


The  radiated  field  can  be  computed  under  the  following  approximation.  The  ground  plane  and  the 
dielectric  are  ignored,  and  a  magnetic  current  filament  is  impressed  in  the  presence  of  the  ground  plane. 
The  magnetic  current  flows  around  the  edge  of  the  patch  and  is  equal  to  the  voltage  between  the  patch  ana 
the  ground  plane.  The  surface  wave  can  be  determined  by  allowing  this  same  magnetic  current  filament  to 


act  in  the  presence  of  a  grounded  dielectric  slab,  although  for  sufficiently  thin  elements,  the  power  carried 
in  the  surface  wave  is  negligible.  The  magnetic  current  is  actually  sampled  at  discrete  points  and  linearly 
interpolated  between  the  sample  points.3  The  radiated  field  is  computed  analytically  and  the  associated 
radiated  power  is  determined  by  numerical  quadrature. 


The  electric  stored  energy  is  computed  by  integrating  the  square  of  the  electric  field.  That  is, 

We  =  I U  J j  |£pdS  *  i U  jjs  |#(r)|2  dS.  (4) 

The  singularities  at  the  load  locations  in  equation  (3)  are  logarithmic;  their  squares  are  integrable.  As 
a  result,  equation  (3)  can  be  used  for  tj£?j  with  little  error.  Since  the  modes  of  the  unloaded  cavity  are 
orthonormal,  the  double  integration  above  reduces  to4 


u)WE  =  51 IlIi  51 

A  r=o <=o  a 


^mn  (<*i')V,mn(**r) 
(** -«.„)*  ' 


Once  the  quality  factor  of  the  loaded  element  is  determined, 


then  its  input  impedance  can  be  computed. 


2.3  Input  Impedance 


Assume  that  one  more  port  is  added  to  the  antenna,  the  ( N  +  1)“  port.  The  circuit  model  as  seen  from  the 
( N  +  1)**  port  for  the  loaded  cavity  is  illustrated  in  Fig.  7.  The  circuit  model  of  the  antenna  is  shown  in 
Fig.  8.  The  -nput  admittance  of  the  lossy  tank  circuit  only  is 


Vr  =  G 
ss  G 


1  +  iQ 


»((£)’- ')] 
V**  (£-')]■ 


where  Q  is  the  quality  factor  of  the  lossy  tank,  and  wo  is  the  resonant  angular  speed  of  the  lossy  tank.  These 
are  given  by 


R  _  2 ujWe 
wqL  P 


where  P  is  the  power  dissipated  in  the  resistor,  and  WE  is  the  time-averaged  energy  stored  in  the  capacitor. 
The  power  P  is  related  to  the  conductance  G  and  the  voltage  V%  across  the  tank  circuit  by 

IK  1 2G=P.  (5) 

The  circuit  quantities  are  related  to  the  antenna  by  choosing  <3  to  be  the  antenna  quality  factor,  and  P  to 
be  the  power  radiated  by  the  loaded-element  resonant  mode  given  by  equation  (3).  The  voltage  Vt  is 

3  See  Appendix  A  for  details  of  this  computation. 

4  The  scries  for  Wg  can  also  be  accelerated,  the  details  of  which  are  given  in  Appendix  A. 
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where  r{  is  the  center  of  the  feed  for  the  loaded  element — the  location  of  the  (N  +  1)**  port.  Note  that  Vt  is 
not  the  feed  voltage  ty  since  the  former  does  not  include  the  effect  of  the  feed  reactance.  The  conductance 
G  can  thus  be  computed  from  equation  (5). 


The  feed  reactance,  A'f ,  can  be  modeled  as  a  linear  process  over  the  band  of  operation.  Let  the  reactance 
corresponding  to  the  lossless  loaded  cavity  ( R  =  0)  be  denoted  by  Xr.  Let  the  total  reactance  be  X(f).  Let 

Xt  =  Xto  +  (/-/o)*n, 
where  /o  corresponds  to  u>o-  Finally,  let 


The  resonant  component  of  the  reactance  can  be  written  as 
_  2xL(fp  +  Af) 
r  1  -  (1  +  A///o)2  ’ 

where  Af  =  f  —  fQ.  For  small  Af,  this  can  be  written  approximately  as 


_ _ rr  fo  _  Tr 

r  Af  T’ 


Then, 


xa*\{X{U)  +  x(f-))  +  ^. 


Similarly, 

An  *  -  [X(U)  -  XU-)  +  2Qw0 L) . 
The  quantity  uioL  =  R/Q  so  that 

*ro»i(A(/+)  +  X(/_))  +  ^, 

and 


XntSfo  [A(A)-A(/-)  +  2/J]. 

What  remains  is  to  determine  the  input  reactance  of  the  ideal  cavity  at  the  upper  and  lower  frequencies  /+ 
and  /_. 


Let  the  indices  of  the  loaded  ports  be  numbered  from  1  to  N  with  corresponding  port  voltages  and 
currents  of  Vi ,  ,  V3,  . . . ,  Vjv  and  /j ,  ,  Jj,  The  (N  +  1)*‘  port  is  the  feed  port  with  feed  current  of 

It  =  1  and  feed  voltage  Vf  =  Z,„  where  Zm  is  the  input  impedance  of  the  loaded  element.  If 


and  I  is  redefined  as 


then, 

V  =  ZI+Zmf/f, 

Vt  =  Z^I  +  Ztt. 

The  impedance  matrix  Z  contains  the  multiport  impedances  between  the  loaded  ports.  The  (1  x  jV)-matruc 
Zmf  contains  the  mutual  impedances  between  the  feed  port  rt  and  the  load  ports  rm.  The  self-impedance 
of  the  feed  port  in  the  unloaded  cavity  is  Zr.  Let  the  diagonal  load  impedance  matrix  be 

ZL  =  Diag  [z£\  Z™ . ZiN) ]  . 

Then 

V  =  -ZLI, 

I  ss  -(Z  +  ZL)-1Zmf,  and 
jX(f)  =  =  Zff-Z|nf(Z  +  ZL)-‘Zmf. 

Having  thus  computed  X(f),  one  can  now  evaluate  the  input  impedance  Z,„  of  the  antenna  using 
Zm  w  Xn  +  Xn(f  -  fo)  +  +  j2Q(f/f0  -  1)] ' 

As  the  results  of  the  next  section  demonstrate,  this  expression  is  reasonably  accurate  over  the  band  of  a 
single  resonant  mode  of  the  loaded  element. 

3  Application  to  Impedance  Matching  and  Results 

In  this  section,  those  elements  of  the  theory  presented  in  the  previous  sections  specific  to  the  application 
of  varying  impedance  are  presented  together  with  corroborating  experimental  results.  The  objective  in  this 
application  is  to  vary  the  input  impedance  while  keeping  the  resonant  frequency  and  pattern  of  the  element 
fixed.  This  objective  can  be  met  by  placing  short-circuits  at  one  or  more  appropriately  chosen  locations 
on  the  patch.  These  short-circuits  can  be  PIN  diodes  thus  allowing  adaptation  of  the  element  to  changing 
conditions  such  as  the  scan  angle  in  a  phased  array.  The  theory  and  results  given  below  assume  that  the 
loads  are  short  circuits.  One  could  probably  also  use  voltage-controlled  capacitors,  varactor  diodes.  Use  of 
varactors  would  probably  simplify  the  biasing  of  the  diodes.  This  possibility  is  currently  under  study. 


Patch 


Figure  9:  The  lociu  of  short-circuit  locations  for  constant  mount  frequency  of  the  single-loaded  element. 


As  seen  in  the  previous  section,  one  can  vary  the  frequency  by  varying  the  location  of  a  short  circuiting 
pin.  It  is  possible  to  vary  the  input  impedance  with  a  single  short  without  altering  the  resonant  frequency. 
The  way  that  the  input  impedance  is  varied  for  a  fixed  feed  location  is  to  move  the  nodal  line  of  the  loaded 
resonant  mode  by  moving  the  shorting  pin  location.  The  closer  the  nodal  line  is  moved  to  the  feed,  the 
lower  the  peak  input  resistance  will  be.  The  trick  is  to  move  the  shorting  pin  in  such  a  way  that  at  least 
the  resonant  frequency  is  unchanged.  This  was  done  using  an  interactive  computer  program  with  graphical 
output  based  on  the  theory  discussed  in  the  previous  section. 

3.1  Single  short 

The  results  of  our  computer  program  indicated  that  if  a  single  short-circuit  were  placed  on  the  locus 
illustrated  in  Fig.  9,  then  the  resonant  frequency  of  the  resulting  loaded  elements  should  be  identical.  Using 
just  a  single  short,  one  can  vary  the  input  impedance  as  is  illustrated  in  Figs.  10-14.  These  figures  contain  the 
computed  magnetic  current  distribution,  the  computed  input  impedance,  and  the  measured  input  impedance 
of  the  loaded  elements.  In  all  cases,  the  patch  dimensions  were  4.00  cm  wide  by  6.00  cm  long  by  0.158  cm 
thick.  The  dielectric  substrate  was  3M  Corporation’s  CuGad  250  glass  reinforced,  double  clad  PTFE.  The 
manufacturer-supplied  nominal  dielectric  constant  was  2.43.  The  feed  point  was  always  at  (1.00, 1.00),  where 
the  coordinates  represent  the  distance  in  cm  from  a  reference  corner  of  the  patch.  The  input  impedance  was 
collected  using  a  computer-controlled,  Hewlett-Packard  8410  network  analyzer  at  the  University  of  Houston. 

It  is  quite  clear  from  Figs.  10-14  that  the  magnetic  current  distributions  around  the  edge  of  the  patch  are 
asymmetrical.  The  pattern  that  results  from  such  a  distribution  will  have  a  significant  cross-polarized  pattern 
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Figure  11:  (4)  Meatured  ud  competed  input  impedance  of  a tingle-loaded  element  with  ihort-circuit  at 
(2.00,1.00).  (b)  Theoretical  magnetic  current  distribution. 


Figure  13:  (a)  Measured  and  computed  input  impadanct  of  a  (ingle- loaded  (lament  with  short-circuit  at 
(2.50,0.91).  (b)  Theoretical  m syndic  current  distribution. 


although  the  .E-plane  pattern  will  be  typical  of  the  (0,  l)-mode  of  a  rectangular  element.  The  theoretical  and 
experimentally  determined  patterns  for  a  patch  loaded  with  a  single  short  at  (2.50,0.91)  are  illustrated  in 
Fig  15.  The  measured  patterns  were  taken  with  the  element  mounted  on  a  1.05m  (8  7-wavelength)  circular 
ground  plane.  The  element  was  driven  at  a  frequency  of  2470  MHz.  Our  pattern  range  is  fairly  new.  it  is 
difficult  to  say  how  much  of  the  cross-polarization  is  due  to  errors  introduced  by  the  range  and  how  much 
is  due  to  the  asymmetry  in  the  element’s  magnetic  current  distribution.  Nevertheless,  the  cross-polarized 
patterns  appear  to  be  at  least  qualitatively  consistent  with  the  predicted  patterns.  Note  that  some  of  the 
discrepancy,  especially  the  scalloping  of  the  pattern  and  the  loss  of  gain  at  the  horizon  is  due  to  the  finite 
size  of  the  ground  plane. 

In  any  case,  it  is  clear  that  the  cross-polarization  arising  when  a  single  short  circuit  is  used  is  significant 
and  probably  unacceptable  in  many  applications.  Because  of  this,  at  least  two  shorts  should  be  used  in 
symmetrical  locations  There  is  another  reason  why  two  or  more  shorts  should  be  used.  It  is  clear  from  the 
measured  and  computed  input  impedances  that  the  impedance  range  obtained  is  not  extremely  broad  for 
a  single-loaded  element.  We  did  not  measure  impedance  levels  for  shorts  any  closer  to  the  feed  than  those 
illustrated  in  Figs  10-14  because  of  the  obstruction  caused  by  mounting  an  SMA  feed  However  we  have 
plotted  the  theoretical  resistance  R  versus  x ,  where  x  is  the  x-coordinate  of  a  single  short-circuit  placed 
along  the  locus  illustrated  in  Fig  9  This  plot  is  contained  in  Fig.  16.  This  result  shou'd  be  compared  with 
the  results  for  the  double-loaded  elements  considered  next. 

3.2  Two  shorts 

With  two  shorts  located  symmetrically  with  respect  to  the  long  center  line  of  the  patch  one  can  eliminate 
the  cross-polarization  in  the  H-plane.  Furthermore,  one  can  more  easily  obtain  a  wider  variation  of  input 
impedance  Figures  17-19  contain  the  theoretical  magnetic  current  distribution,  and  the  measured  and 
computed  input  impedances  of  elements  loaded  by  symmetrically  located  pairs  of  short-circuits.  The  locus 
of  points  on  which  these  pairs  of  short-circuits  should  lie  is  illustrated  in  Fig  20. 

The  theoretical  and  measured  patterns  are  given  in  Fig.  21.  The  measured  cross-polarization  is  dis¬ 
appointingly  high  in  this  case.  As  mentioned  before,  this  may  be  due  to  errors  introduced  by  the  range. 
Nevertheless,  the  cross-polarization  is  definitely  reduced  from  what  it  was  in  the  case  of  the  single-loaded 
element. 

Finally,  it  is  clear  from  Figs  17-19  and  Fig.  22  that  the  impedance  variation  is  greater  for  the  double- 
loaded  element  than  for  the  single-loaded  element.  In  short,  it  appears  that  the  performance  yielded  by  a 
single  short  falls  short  of  that  obtainable  using  two  shorts. 


4  Conclusions  and  Future  Research 

There  are  several  conclusions  that  can  be  drawn  from  the  theory  and  results  presented  in  this  report. 

1.  It  is  possible  to  change  the  input  impedance  of  a  microstrip  element  over  a  wide  range  without  affecting 
its  resonant  frequency  or  co-polar  pattern  by  moving  short-circuits  from  one  point  to  another. 


Figure  15:  (a)  Measured  and  computed  input  co- polar  and  croes-polar  patterns  in  the  JF>plane  for  a 
short-circuit  at  (2.50,0.91).  (b)  Measured  and  theoretical  patterns  in  IT-plane. _ _ 


0>) 

Tifort  18:  (a)  Measured  and  computed  input  impedance  of  a  double-loaded  element  with  short-circuits 
at  (2.50, 1.21)  end  (2.50,2.79).  (b)  Theoretical  magnetic  current  distribution. 


2.  The  theory  that  predicts  this  result  has  been  corroborated  by  experiment.  Specifically,  the  experiment 
shows  that  the  resonant  frequency  of  the  mode  is  left  relatively  unchanged  by  the  appropriate  placement 
of  short-circuits.  Small  deviations  between  the  predicted  and  observed  resonant  frequency  do  exist. 
They  are  due  partly  to  errors  in  fabrication  of  the  element  and  partly  due  to  the  errors  introduced  by  the 
cavity  model  upon  which  the  theory  is  based.  The  experiment  also  shows  that  the  agreement  between 
the  predicted  and  observed  input  impedance  is  excellent  in  most  cases.  Finally,  there  was  qualitative 
agreement  between  theoretical  patterns  and  measured  patterns  although  the  measured  cross-polarized 
components  in  the  double-loaded  element  were  much  higher  than  expected  from  the  theory—- possibly 
a  problem  with  our  range  and  ground  plane  fixture. 

3.  From  the  point  of  view  of  producing  lower  cross- polarization  and  from  the  point  of  view  of  more  easily 
varying  the  impedance  level  over  a  wide  range,  two  short-circuits  are  better  than  one  short-circuit. 
Each  short-circuit  pair  should  be  biased  together  with  the  same  biasing  voltage. 

Some  important  items  are  left  for  future  research.  These  include  the  following. 

1.  A  simulation  of  loaded  elements  in  a  phased  array  needs  to  be  earned  out  to  determine  the  basic 
feasibility  of  the  ultimate  application.  This,  we  propose,  as  the  major  deliverable  of  a  future  research 
contract. 

2.  Because  each  diode  or  diode  pair  needs  a  biasing  line,  the  number  of  biasing  lines  could  become  quite 
large  if  many  levels  of  adjustment  in  impedance  are  provided.  For  this  reason,  we  are  currently  studying 
the  feasibility  of  using  voltage-controlled  capacitors  to  change  the  impedance  level. 


A  Additional  Theoretical  Results 


This  appendix  contains  additional  theoretical  results  that  are  implemented  in  the  computer  program  used 
to  generate  the  data  presented  in  Section  3.  There  are  three  topics  within  this  appendix: 

1.  The  efficient  computation  of  the  cavity  Green’s  function, 

2.  The  computation  of  the  far-held  pattern,  and 

3.  The  computation  of  the  energy  stored  in  the  cavity. 

A.l  The  efficient  computation  of  the  cavity  Green’s  function 

The  expression  for  the  Green’s  function  in  equation  (1)  is  very  slowly  converging  and  not  particularly  useful, 
especially  for  patches  of  other  than  rectangular  shape.  It  can  be  accelerated  in  the  following  way.  First 
partition  the  set  A  of  mode  indices  into  the  union  of  three  sets, 

A  =  {(0,0)}  Ul*  UA(. 

The  first  set,  consisting  of  a  single  element,  represents  the  indices  of  the  “D.C.  mode.”  This  is  the  mode  of 
the  unloaded  cavity  that  resonates  at  zero  frequency — it  is  just  a  constant.  Since  the  modes  are  assumed  to 
be  normalized,  this  constant  is 

where  S  represents  the  area  of  the  patch. s  The  third  set  is  the  set  of  modes  whose  resonant  wavenumbers 
are  greater  than  or  equal  to  (fcmaL<: 

A(  =  {(m,  n)|t  mn  mmx  }  • 

The  parameter  fcmu  is  the 

L  _  2x/mu^/£r/ir 

] 

C 

where  c  is  the  speed  of  light  in  a  vacuum,  and  fmmx  is  the  maximum  frequency  of  interest  for  the  antenna — a 
number  that  can  be  specified  a  prton  for  any  analysis.  The  parameter  {  is  some  real  constant  greater  than 
one.  The  second  set  is  the  complement  of  A(  with  respect  to  the  set  Aq, 

A(  =  Ao  -  A{ , 

where 

>10  =  A  -  {(0,0)} 

is  the  set  of  all  mode  index  pairs  except  for  that  corresponding  to  the  D.C.  mode. 

Suppose  that  |<|  <  l/{.  Then  one  can  approximate  the  fraction  1/(1  -  £J)  by 


-  _  *ao+  a iCJ  +  a2C*  +  ■  •  •  +  at_1CJ(i_1). 

Then  for  mode  indices  in  A(,  one  can  approximate  the  fraction  l/(Jfc2  -  fc*,B)  by 


k7  -  k 7 

K  fcmn 


1  k7  k7L~ 7 

~  -aoii — a'u - 

“mn  'mn  ftmn 


Thus,  the  sum  in  equation  (1)  can  be  partitioned  and  approximated  as 


G(r|r')  »  ju>p 


......  J  ^  ,  V'  *l>mn(,r)ll>mn(r')  .  V"  V’mn(r•)V’mn(»■,) 

G(rir  )  «  ^  +  2_  LT_  l; - “o  2- - 15 - 

l4  A,  m" 

l2  'i>mn(r)tl>mn(r ')  ,  JL_2  ^  tl>m„(r)^mn(r')  j 

-“1*  2- - M - 2- - pE -  f 

A(  mn  At  mn  ) 

This  sum  can  be  recast  in  a  more  convenient  form  for  acceleration  by  summing  over  Ao  instead  of  A^  The 
result  is 

G(»V)  *  J'wp  <  in  H  '  +aoJ-  +  ai]T~-b--  -bajr_1i:5I-1  ^mn(r)V-mn(r') 

y  ^  —  L*  *mn  *mn  *mn  Kmn  J 

_ffl02- - p - a»*  2- - p - 

A,  mn  A,  mn 

- «.,» «-»£i=S££d2\. 

a,  *mn  J 

There  are  compromises  to  be  made  in  the  choice  of  the  parameter  (.  If  £  is  chosen  too  near  1,  then  the 
required  number  of  terms  L  necessary  to  approximate  the  fraction  1/(1  —  (7)  becomes  very  large.  On  the 
other  hand,  if  (  is  chosen  to  be  much  greater  than  1,  L  may  be  small,  but  the  number  of  elements  in  the  set 
A(  becomes  very  large.  Either  of  these  conditions  is  undesirable.  A  reasonable  choice  of  £  is  £  =  2  If  L  is 
then  chosen  to  be  2,  and  the  coefficients  ao  and  ai  are  chosen  to  be  0.989  and  1.331,  respectively,  the  error 
in  the  approximation  is  about  one  percent.  Thus,  the  Green’s  function  can  be  approximated  as 

G(»V)  =  L2  _*  L2  +aon~-Min-  V’mn(r)V>mn(r')  (6) 

y^  -  l*  “  *mn  *mn  Kmn. 

+  aoGo(r|r')  +  fcJ<JiGi(r|r')^  , 


+  aoGo(r'|^,)  +  ^aiG^i'lr' 


where 


Go(r|r')  =  -£^ 


& 


GiMO  =  -E 

Ao 


^mnW^mnv  ) 


The  decomposition  of  the  Green’s  function  expressed  in  equation  (6)  represents  an  expansion  of  G  in 
terms  of  its  first  few  lowest-order  modes  plus  the  contribution  due  to  ail  higher-order  modes,  a  quasistatic 
expansion.  Once  the  quasistatic  terms  Go(r)r')  and  Gi(r|r')  are  computed,  then  the  Green’s  function  can 
be  evaluated  very  efficiently  for  different  frequencies — the  most  difficult  part  will  have  been  done  once  and 
for  all  for  a  given  pair  of  source  and  observation  points  r  and  r‘. 

The  series  for  Go(r|r')  and  Gi(r|r')  are  still  very  slowly  converging,  as  slowly  converging  as  the  series 
for  G.  These  quasistatic  Green’s  functions  satisfy  the  boundary-value  problem, 

VJG0(r|r')  =  6(r-r')-|,  r  €  S, 

VG0(r|r')-n  =  0,  r  €  8S\ 

VlGi(r|r')  =  -G0(r|r'),  r  €  S, 

VGi(r|r')-n  =  0,  r  6  dS. 

One  can  obtain  Go  and  G\  by  evaluating 

G"  *  <7> 

G'  “  i&{ap(G-F5)}-  W 

By  expanding  G  in  a  Laurant  series  in  k 3  about  k3  =  0,  one  can  evaluate  these  limits.  To  this  end,  the  patch 
is  first  described  in  terms  of  a  coordinated  system  appropriate  to  the  patch  shape. 

Let  u  and  v  be  the  coordinate  variables  of  an  orthogonal  coordinate  system  derived  through  a  conformal 
mapping  of  the  rectangular  coordinate  system.  This  mapping  is  chosen  so  that  it  maps  the  patch  into  a 
rectangular  patch  in  the  (u,  v)-plane.  Of  course,  if  the  patch  is  already  rectangular,  then  u  =  x,  v  =  y.  For 
a  circular  patch,  a  circular  sector  patch,  an  annular  patch,  or  an  annular  sector  patch,  the  transformation 
used  is 

x  —  oe®  cosu, 
y  =  ae*  sinu, 

where  a  is  the  outer  radius  of  the  patch.  The  transformation  need  not  be  analytic,  but  it  must  satisfy 
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If  these  equations  are  satisfied,  then  the  boundary  value  problem  for  G  is  transformed  into 
[V2  +  fc2/»(u,  t>)2]G(r|r')  =  ju^6(r  -  r'), 

where  h(u,  v)  is  the  common  scale  factor  of  the  u  and  v  variables.  If  the  patch  is  separable,  then 
h(u,v)J  =  /(u)  +  g(v). 

In  this  case,  G  can  be  expressed  as 


where 


(r|  1  - KwWxT) - • 

Nm=  r*Um(u)7du, 

JUm> 

w  (V-,V+)  =  Wm  =  V-(v+)V+'(v+)  -  V-'(v+)V+(v+l 


V, 

V  <  v' 

V1, 

V'  <  V 

V, 

v'  <  V 

v\ 

V  <  v1 

The  prime  “/  ”  on  dependent  variables  indicates  differentiation  with  respect  to  argument.  The  primes  on 
u  and  t  indicate  the  source  point  location  («',  t/)  in  the  (u,  v)-coordinate  system.  The  Um’ a  and  Vm’s  satisfy 
the  separated  Helmholtz  equations, 

U^  +  {aim  +  kif(u))Um  =  0, 

C +  (-«2. +  *’»(«))  ^m=  0, 

where  the  “±”  superscript  has  been  suppressed  on  the  “V  for  convenience. 

The  a  is  the  separation  constant  determined  by  applying  the  pair  of  homogeneous  boundary  conditions 
that  Um(u)  satisfies  at  u  =  u± ,  where  u_  is  the  left  boundary  of  the  transformed  patch  and  u+  is  its  right 
boundary.  Similarly,  is  the  lower  boundary  of  the  transformed  patch  and  v+  is  its  upper  boundary.  These 
boundary  conditions  depend  on  the  geometry  of  the  patch.  They  are  either  the  physical  boundary  condition 
that  the  derivative  Um(u)  must  vanish  at  one  or  both  of  the  limits  of  u,  or  that  Um(u)  must  satisfy  conditions 
that  will  make  the  field  and  its  derivatives  continuous  at  the  limits  of  u.  Without  loss  of  generality,  one 
can  always  assume  that  V£'(v+)  =  0.  The  condition  that  (v)  satisfies  can  also  be  the  physical  boundary 
condition  that  its  derivative  vanishes  at  v  =  ,  or  that  the  field  be  continuous  at  t>  =  v_ . 

These  parameters  can  be  expanded  in  a  power  series  in  fc2  as 

Um  =  Uom  +  k7UUn+k*U3m  +  ---, 


A 


V'm  =  V'0m  +  k7Vim  +  k*V2m  +  ■■■, 

ol  =  <‘L  +  ^L  +  ^L  +  -, 

Am  —  Nom  +  k7Nim  +  k*  N2m  +  •  •  • , 

Wm  =  W'om  +  k7Wlm  +  k*W2m  +  ■■■. 

The  6rom’s  and  the  V^’s  satisfy  the  sequence  of  differential  equations, 

n  —  1 

^nm  “1"  ^Om C"nm  —  ^  ^(n  —  — 


^ nm  aOm^im  —  ^  ^  ^(o  — }(’,)^(n-l),mi 
<s0 

where  again  the  superscript  “±”  on  the  V’s  has  been  dropped  for  convenience.  In  addition  to  applying  the 
boundary  conditions  discussed  above,  the  following  conditions  are  also  applied  to  obtain  the  solution  in  the 
most  convenient  form: 

Aim  —  Ajm  es  •  •  =  0, 

=  0,  n  >  0, 

*£.(«+)  =  1, 

KTm(v- 1  =  0,  n  >  0, 

Vom(^-)  =  l- 

These  conditions  can  be  applied  since  any  of  the  multiplication  of  the  U  and  V  functions  by  any  arbitrary 
function  of  k1  leaves  G(rjr')  unchanged,  but  does  change  the  coefficients  of  the  series  in  k7  listed  above — 
i.e.,  the  expansions  given  above  are  no<  unique.  This  choice  simplifies  the  coefficients  of  the  expansion  of  the 
Wronskian,  they  become 

W'nm  =  -V^{V+). 

Using  the  solutions  of  these  equations  in  equations  (7),  (8),  and  (9)  yields  the  expressions 


G<(*-|r')  =  rl0  +  ^2  r,m. 


The  r  s  are 


-H-i-W* 

=  |  |n20  -  n^W20  -  ^30  +  (^K>)  |  , 


'  ^ *.*/*-" 


where  Au  =  u+  —  ti_.  The  n’s  are  given  by 

nom  =  C'om(u)tr0m(tt,)V^(w<)V'c+(t»>), 

nim  =  Uim(u)U0m(u,)V0-m(v<)V+m(v>)  +  C/om(u)y1m(u')V0-m(v<)V0+ (v>) 

+  ^0m(«i)C/om(u')Vfm(«)<)V0+m(t,>)  +  t/om(«)C/om(«')K>m(«<)nm(t'>). 

n20  =  Cr2o(«)  +  ^2o(«')  +  ''Vo  (v< )  +  mo(v>) 

+  U\o{u)Uio{u')  +  Vio(w<)^io(v>) 

+  Vio(tt)Vio(»<)  +  t/io(«)V+(t,>)  +  UMv'WfoM  +  <7io(u,)V+(v>). 

The  summand  of  Go  always  contains  the  term 


where  a0m  =  mi/,  z  —  v>  —  v<  +  j(u  —  u').  This  term  can  be  summed  in  closed  form  to  yield 

Go(r|r')  =  |i/t'>  2"t)<  -  In  \vz\  -  In  Sm^yy—  j  +  remaunder  terms. 

For  source  points  (u\  v')  not  on  the  edge  of  the  patch,  the  remaunder  terms  of  Go  in  the  preceeding  expression 
are  non-singular  as  (u,  tz)  —  (u\  tz')  The  logarithmic  singularity  of  the  Green’s  function  is  explicitly  extracted 
by  this  method.  The  remainder  terms  of  Go  and  of  Gi  involve  sums  over  m  of  terms  of  the  form 

(mi/)' 

where  r,  is  a  complex  number  whose  real  and  imaginary  pauts  are  linear  functions  of  u,  ti\  t/,  and  iz',  and 
t  =  1,2,3,  .. 

The  senes  associated  with  these  terms  cam  be  evaluated  efficiently  using  the  functions 
Fn(r)  =  K  **1,2....,  £{r}>0. 

The  first  function  Fx(z)  can  be  summed  in  terms  of  elementary  functions  and  is  given  by 
*i(*)=  ^-ln{2|sinh(z/2)|}. 

The  second  and  third  functions  Fz(z)  and  Fz(z)  can  be  computed  for  z  in  various  regions  of  the  complex 
r-plane  as  indicated  in  Figs.  23-24. The  formulas  are  extended  by  the  following  identities  to  regions  of  the 
complex  r-plane  not  covered  in  these  two  figures: 


34 


Imaginary  axis 


Figure  24:  Formal*!  for  Fj(»)  for  different  region*  of  the  complex  j-plene. 


/^(z)  =  F„(zmod2T) 
Fn(z)  =  Fn(z’). 


This  method  has  been  used  to  obtain  explicit  formulas  for  the  several  differently  shaped  patches.  The 
formulas  for  rectangular  patches  are 


*o  =  -[(y>  -  y<)  +  i(*  -  *')], 

a 


z\  =  —  [(y>  +  y<)  +  j(*  -*')]. 


*2  =  —{[2*  —  (y>  +  y<)]  +  j(* -*')}. 

a 


=  -l(y>  -  y<)  +  j{*  +  *')]. 

a 


*<  =  -[(y>  +y<)  +  ;(*  +  *')], 


z5  =  —{[26  —  (y>  +  y<)]  +  j(* +  *')}. 


r°°  =  S  ^ +  (*  —  . 

no  =  ^{^<  +  (6-v>)2]  +  v2<(i-y>)2-^V<  +  (fc-y>)2]  +  ^2}. 

=  -»{*<},  <  =  0,1,2, 

<7-4  =  —~[26  —  (y>  —  y< )], 
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coth(mT6/a)  fWylyO  -  ftim(yly'), 
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Gx(r\r')  =  r10 


-£{£>  [«w*>  - 1  (i  -  •»)  «”•■]  -  (==) « (=?) 


These  series  typically  do  not  require  more  than  two  or  three  terms  for  adequate  convergence. 


A. 2  Computation  of  the  Far-Field  Pattern 

Given  a  patch  with  a  set  of  L  sample  points  (ti/,  «/)  along  its  edge  (with  (ti£,  Vi)  —  («i,  t>0),  where  (u,  v)  is 
a  point  in  the  curvilinear  coordinate  system  used  to  describe  the  patch,  one  can  compute  the  approximate 
far-field  in  the  following  way.  First,  sample  the  loaded  element  mode  <f(u,  v)  at  each  of  the  sample  points 
thus  forming  the  set 

Then  the  far-field  components  of  the  electric  field  are  given  by  the  following  formulas. 

Ei  =  jk0t(FT  sin  <t>  —  Fy  cos  4>)——c~*ktr, 


where 


Ei  =  jk0t(Fx  cos  <t>  +  Fyain4)cos6-^—e~>kar, 


F  =  jo •  Ar/^  +  j  ji  •  Ar<)]  A*/e>t,,r  r', 

=  £[*(wr+i.«'/+i) +  *(«/,«/)], 

=  *(u<+1,vr+l)  -  ♦(w,  vt), 

rt  =  */f  +  y/y,  F/  =  i(r/+i  +  r<),  Ar/  =  rf+l  -  rt, 
m  t  =  u(xt,iu),  t>/ =  w(x/,y/),  As/  =  |Ar/|,  f  /  = 

415/ 

f  =  x  cos  ^  sin  +  y  sin  ^  sin  5  +  zcosfl, 


Jo(i)  = 


ii(*)  = 


sin  z  cosz 
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A.3  The  Computation  of  the  Energy  Stored  in  the  Cavity 

This  section  deals  with  the  efficient  computation  of  the  electric  stored  energy  within  a  microstrip  antenna 
loaded  with  N  lumped,  linear  loads.  The  expression  for  the  loaded  mode  distribution  is  given  in  equation 
(3)  in  Section  2.  The  actual  load  currents,  which  are  distributed  over  a  finite  area,  are  assumed  to  be  fila¬ 
mentary  for  the  purpose  of  computing  energy.  This  can  be  done  because  while  the  field  itself  arising  under 
this  assumption  has  log  singularity,  the  square  of  this  field  over  the  volume  of  the  cavity  is  integrable.  The 
associated  stored  energy  times  u>  is  given  by  equation  (4)  Section  2.  This  series  converges  slowly  enough 
that  some  moderate  acceleration  is  useful.  The  series  is  summed,  approximately,  by  extracting  a  single 
quasi-static  term.  This  quasi-static  term  is  identical  to  the  1 '‘-order  quasi-static  term  of  the  expansion  of 
the  Green’s  function  itself. 

The  ratio 

[t,-V  (>») 

V*  *mn)  'mn  *mn 

for  (m,  n)  €  The  set  A(  is  the  same  as  that  defined  in  Section  A.l.  For  the  sake  of  numerical  efficiency, 
the  same  value  of  £  should  be  used  in  this  computation  as  that  used  for  computing  the  field  in  equation  (6) 
of  Section  A.l.  Using  the  approximation  in  (10),  one  arrives  at 


.  N  N 

WB  =  Jetn3rk3orl0YtY,ItI< 


where 


■  I  s* + ?  -ku*  ■  sL  ■ 61  ^»(r<)^n(r/) 
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-  AoGi(ri|r*)  —  6ik2Gj(ri|r*)|  , 

G,(r|r')  =  -J2  </,mn(r)6t/'m"(r,). 
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The  computation  of  Gi(r|r')  is  discussed  in  Section  A.l.  The  computation  of  G3(r|r')  could  be  handled 
in  the  same  way  as  the  computation  of  the  lower-order  quasistatic  terms,  but  this  is  too  tedious  a  computa¬ 
tion.  Instead,  Gj(r|r')  is  summed  directly. 

The  coefficients  6o  and  6;  are  the  solution  of  the  minimax  problem, 

min  <  max  |»7(x)|  >  , 


A  plot  of  rj(x)  versus  x  will  show  that  the  minimax  solution  will  satisfy 

m  =  *?(i/£2) 

and 

»?(0)  =  ~  »?(*mix). 

where  xmtx  satisfies 

d  .  .  2  . 

sr**' =  ‘ 

The  solution  to  this  is 


The  solution  to  the  minimax  problem  is 


(1  -  1/^)2  -  1 


For  i  —  2,  these  parameters  are 
60  =  0.958 

and 

6t  =  3.111. 

This  approximation  yields  a  maximum  error  of  about  4%. 

The  next  problem  is  to  find  Gj(r  |r').  WhVt  is  the  error  if  Gj(r|r')  were  only  summed  over  Aj?  If  £  =  2, 
then  a  typical  relative  error  level  would  be  no  worse  than 

2-“  «  0.06. 

This  is  because  kmtx  is  typically  greater  than  or  equal  to  the  smallest  non-zero  kmn.  Contributions  to  the 
sum  outside  of  A;  are  at  least  2~e  times  as  small  as  the  contribution  due  to  the  smallest  kmn.  In  the 
worst  case,  all  contributions  outside  of  A(  add  in  phase  (when  the  source  and  observation  points  coincide) 
in  which  case  the  sum  can  be  approximated  by  a  continuous  integral  over  wavenumber.  This  accounts  for 
the  exponent  of  4  in  the  previous  equation.  This  approximation  of  truncating  Ga(r|r')  to  just  a  sum  over 
A(  should  be  adequate.  Under  this  approximation,  the  expression  for  We  becomes 
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Command,  Control,  Communications  and  Intelligence 
( C 3 1 )  activities .  Technical  and  engineering 
support  within  areas  competence  is  provided  to 
ESV  Program  Offices  (PCs)  and  other  ESV  elements 
to  perform  elective  acquisition  0($  C3I  systems. 

The  areas  o  6  technical  competence  include 
communications ,  command  and  control,  battle 
management,  information  processing,  surveillance, 
.sensors,  intelligence  data  collection  and  handling, 
solid  state  sciences ,  electromagnetics ,  and 
propagation,  and  electronic,  maintainability , 
and  compatibility . 


